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Abstract. In the dissipative quantum dynamics of a mesoscopic aggregate of excited two 
level systems (atoms) coupled to a single resonance mode of a cavity, two physical phenomena 
associated with superradiance appear. A pronounced emission peak on short time scales is 
related to the known superradiant burst of excited atoms cooperatively radiating into free 
space. It is followed by relaxation to a stationary state of the composite system such that 
for sufficiently large atom-field coupling a strongly correlated state emerges. The crossover 
to this state can be interpreted as the precursor of the transition from normal to superradiant 
phase in the thermodynamic limit of the Dicke Hamiltonian. Motivated by recent experimental 
activities, these features are investigated in detail for a mesoscopic number of atoms and a 
cavity embedded in a dissipative medium described by a damped Dicke model. We identify 
observables and characteristics of the quantum dynamics on shorter time scales which allow 
to clearly distinguish weakly correlated from strongly correlated Dicke physics. 
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The Dicke model, originally developed to describe a large number N of two-level atoms 
interacting with a single-mode radiation field [1, 2], has regained substantial interest in the 
last decade [3-6]. Recent experimental realizations of Dicke physics have utilized various 
degrees of freedom of cold atom clouds, see e.g. [4, 5, 7], and implementations with solid 
state devices have also been discussed [8]. 

In this context, theory has contributed to a variety of aspects of the model with most of the 
work basically focusing on two issues, namely, the dissipative dynamics in the semiclassical 
limit of a large collective spin consisting of the individual atomic levels, and properties at and 
close to the equilibrium state in the thermodynamic limit of ./V —> °° and vanishing coupling 
to a dissipative medium. Both regimes display characteristic phenomena associated with 
the term superradiance: on a transient time scale, the cooperativity of the atoms may lead 
to superradiant light emission (burst) while thermodynamically a phase transition occurs 
when the atom-field coupling exceeds a certain threshold. The latter threshold captures 
the changeover from a regime where energy eigenstates are product states of the atoms 
and the field mode, to a regime where collective behavior associated with strong atom-field 
correlations appears (superradiant phase). 

What has gained much less attention, is first the relation between these two facets, the 
superradiant burst and the phase transition, with respect to atom-field correlations. This issue 
is not only relevant experimentally [5], but also for the analysis of non-equilibrium phase 
transitions [9]. Second, the impact of dissipation and finite size effects on the emergence 
of a superradiant phase have not been fully understood yet, despite their important role for 
mesoscopic aggregates. The purpose of this work is to contribute towards closing this gap. In 
particular, we will focus on the non-equilibrium dynamics of the Dicke model starting initially 
far from equilibrium for a moderate number of atoms (mesoscopic regime). 

In fact, experimental progress for tailored superconducting circuits allows to integrate 
multiple artificial atoms in the form of Cooper pair boxes into micro-cavities. To maintain 
sufficiently long coherence times and tunability, the number of two-level systems in these 
devices must be kept moderate and far below the number of atoms in a cold atomic cloud. 
Furthermore, environmental degrees of freedom in the form of heat baths are always present 
in these solid state structures. Nevertheless, one may assume that signatures of Dicke physics 
are observable [6], at least in certain ranges of time and parameter space. Another line of 
research combines atomic ensembles with superconducting cavities to benefit from the long 
coherence times of the former and the fast tunability of the latter entities [10]. In this case, 
dissipation on the solid state side, i.e. the finite photon lifetime in the micro-cavity, is the 
dominant mechanism and must be taken into account. 

These new composite systems in turn trigger the more fundamental question, namely, to 
what extent the Dicke model is realized in these aggregates at all. For example, the atom- 
field interaction may include higher order terms in the field operator not captured in the 
model [8] or couplings to residual degrees of freedom may be present (e.g. trapping loss 
of atoms, charge-background fluctuations in solid state systems, etc.). The Dicke model 
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may be then justified at least on transient time scales which are short compared to typical 
relaxation times of the full composite system. Consequently, it may be extremely difficult, 
if not impossible, to monitor Dicke dynamics on long time scales and towards relaxation to 
a collective thermodynamic ground state. Here, we thus also address the question whether 
within a mesoscopic Dicke aggregate, signatures of the true phase transition are encapsulated 
in superradiant-type phenomena on transient time scales and in the presence of dissipation. 

The paper is organized as follows. We start in section 2 with a brief discussion of the 
two characteristic phenomena, namely, the superradiant burst and the superradiant phase and 
then introduce the dissipative Dicke model explored in the sequel. In section 3, the dynamics 
of the photon population in the field mode (cavity) is studied in detail within the full model 
which is contrasted with its often used form based on a rotating wave approximation. Scaling 
properties of the initial radiation peak are analyzed. The issue of how to understand and 
monitor experimentally possible collective dynamics of atoms and cavity is addressed in 
section 4. 

2. Two facets of superradiance 

To set the stage, we briefly recap in this section the characteristics of the two well-known 
facets of superradiance and discuss how we envision their interplay in the proposed scenario. 

2.1. Light emission into free space: The superradiant burst 

A single excited atom in free space will radiate its energy isotropically in an exponential decay 
process. If a large number of excited atoms are brought into close proximity their decay is 
not described by adding the light emission rates from all the atoms independently. Instead 
a cooperative emission process of all atoms is established, which results in a strong, short 
anisotropic burst of light. This superradiant free-space burst, predicted by Dicke in 1954 [1] 
and first experimentally observed by Skribanovitz [11], has since been extensively studied and 
its characteristic features are well understood (see e.g. [2] for a review). 

For the purpose of this paper, we will consider as defining signature of superradiance the 
crossover from an exponential decay for a single atom to a peak in the time-dependence of the 
emission rate. This peak gradually becomes more pronounced when the effective number of 
atoms, N, cooperating in the emission process increases (roughly speaking, N is the number 
of atoms within a volume, Xy, where A v is the wavelength of the emitted light). Due to 
the cooperative nature of the light emission, the peak’s height grows as N 2 while its temporal 
width shrinks as A -1 , since the integrated emission is proportional to the total initial excitation 
energy and, hence, to N. 

The simplest theoretical description of such burst features relies on a rate model. 
The radiation of N excited indistinguishable two-level atoms is calculated by determining 
decay rates between the N + 1 levels of a spin-A/2 system (which results from adding the 
(pseudo)spins of two-level atoms, see below) from the corresponding dipole matrix elements 
[ 2 ]. 



Non-equilibrium dynamics of the Dicke model 


5 



(b) 



Figure 1 . Visualizations of the potential landscape of the Dicke Hamiltonian. 

(a) The Holstein-Primakoff approach maps cavity and spin degree of freedom to (displaced) 
bosonic operators (cf. (2) and (3)). Above a critical coupling strength two equivalent potentials 
V ± quadratic in the corresponding (dimensionless) position variables X a / b are found, which 
describe the system close to two mean field solutions (X a , X b ) = (±v^a,Ty^jS), see (5). 

(b) In a J* eigenbasis the Dicke Hamiltonian can be visualized as shifted parabolas for the 
various M x eigenvalues, that are coupled by a J z term, see (6). 


2.2. Thermodynamics: The superradiant phase transition 


2.2.1. The Dicke Hamiltonian A collection of N identical atoms in a single mode cavity can 
be described by the Dicke Hamiltonian [1-3], 


^Dicke — G)(f a + £j' + — y= + aj (J+ -\- J ) , 


( 1 ) 


where the atoms are modeled as two-level systems with identical level spacing e for which 
collective ("large”) spin operators are introduced, J z = ^ Ejli <7?, , h = 1. 

The collective spin couples to a common, single electromagnetic mode of frequency CO in the 
cavity described by bosonic operators a, a ' with effective coupling strength X / s/N. The latter 
one is conveniently introduced when the density of atoms per unit volume is fixed since then 
the bare coupling A between any single atom and the atomic mode is effectively reduced with 
a growing number of atoms. If instead in a given experimental situation the actual density of 
atoms can be increased, A grows accordingly. 

The above modeling (see [3] for a recent review) relies on the assumption that all atoms 
couple to the mode with equal strength (cf. e.g., [12,13] for inhomogeneous effects) and that a 
loss of atoms from the collective state (e.g. due to trapping losses or external noise sources, see 
[14]) can be neglected on the timescales of interest. Accordingly, the collective spin operators 
in (1) couple only states within one J manifold of the so-called Dicke states | J,M). For 
example, these are states with maximum angular momentum J = N /2 when the atomic system 
is initialized in the atomic ground state |7 = N/2,M = —J) or the maximally inverted state 
|/ = N/2,M = +/). The validity of the Dicke Hamiltonian for various experimental scenarios 
with (artificial) atoms has been widely discussed recently [8, 15], Within a rotating wave 
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approximation (RWA), where only the coupling terms PJ~ + aJ + are kept, the Hamiltonian 
conserves in addition to the spin magnitude J also the total number of photonic and atomic 
excitations N ex = d'a+J 7 +J. The full (non-RWA) Hamiltonian (1) only conserves the parity 
n = exp (\nN ex ). 

2.2.2. Phase Transition Considering the ground state of the Dicke Hamiltonian, it is 
immediately obvious, that there is a competition between the coupling term and the terms 
describing the energy of the uncoupled cavity and spin system, respectively. This competition 
depends on the strength of the coupling A compared to the excitation energy for a cavity 
photon ft) or a single spin excitation e. For weak coupling, the ground state of the composite 
system lies close to a product state of the individual ground states of isolated cavity and 
isolated spin, i.e. |GS^ =0 ) := \n = 0, M = —A/2). Of course, this state is a common 
eigenstate of the cavity number operator h = d'd and J z with h GS^ = „)=Oand/'|GS A=0 ) = 
-(Af/2)|GS A . 0 ). 

In the opposite regime of strong coupling, one has a finite expectation value for the 
operator (df + a) ( J + +J~) associated with a strongly correlated ground state with non-zero 
photon occupation and a finite number of spin excitations. In fact, in the thermodynamic 
limit N —) oo, the Dicke Hamiltonian was shown [16,17] to exhibit a phase transition between 
a superradiant phase with macroscopic occupations in the field and the atoms and a normal 
phase without excitations at zero temperature. This transition occurs at a critical coupling 
strength Ac = \fofe /2 [17] and at A C RWA = \fcoe for the RWA version of the Dicke Hamiltonian 
[16]. 

The most intuitive picture for the emergence of this transition and the nature of the 
respective ground states is provided by the Holstein-Primakoff (HP) approach [18]. The HP 
treatment of the Dicke Hamiltonian essentially consists of (see [18] for details): 

(i) Introducing bosonic operators b, P to replace the large spin operator in a HP transformation 

j z = Pb - /, /+ = P VzJ^Pb, J- = y/lJ^Pbb . ( 2 ) 

(ii) Introducing displaced operators 

a — y c zb b — y d -j- \/~f , (3) 

with real-valued parameters a,f5 to be determined self-consistently. In the superradiant phase 
they are assumed to be of order ff(N), corresponding to a macroscopic number of excitations. 

(iii) Performing the HP approximation, a lowest order expansion in the scaled number of 
excitations Pb/N assuming only small fluctuations around the mean field solutions. This 
leads to a Hamiltonian bilinear in the bosonic operators c, d , while linear terms are eliminated 
by the self-consistent choice for the displacements (Bogoliubov transformation). 

It then turns out that for A < A c there is only a trivial solution a = 0 = /3, while for 
A > A c the displacements 

u-tAcM) 2 ] 


(4) 
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indeed correspond to macroscopic expectation values 

= (jhse = P _ . = _yt 

J (axf ’ J J V- ’ 

with (-)sp denoting expectation values in the superradiant phase. 

The respective Hamiltonians above and below the transition can also be formulated 
in terms of two sets of conjugate operators X a /^, P a /i :) , i.e. positions and momenta for the 
respective operators of cavity and atomic degrees of freedom. The position degrees of freedom 
X a /i, are coupled via bilinear potential terms which provides a very intuitive understanding of 
the emergence of a collective ground state. Namely, the principal axes are rotated with respect 
to the original cavity and atomic direction such that the potential surface becomes unstable in 
one direction when A exceeds the critical value Ac. Correspondingly, the HP approximation 
based on the assumption that fluctuations remain small around the fixed points a = /3 = 0 
fails. It only applies in the normal phase A < Ac sufficiently away from the phase transition, 
where then the Dicke model is well described by a two-dimensional harmonic oscillator 
potential centered around X a = X/, = 0. Instead, in the superradiant phase A > Ac the HP 
approach leads to two quadratic potentials V ± centered around either one of the solutions, 
(Xa, Xb) = (±V2 oc,tVW) [see figure 1(a)]. 

In this latter regime, eigenstates of the model which respect the parity symmetry of 
the original Dicke Hamiltonian, are constructed from (anti)symmetric superpositions of the 
individual states of H ± . In the thermodynamic limit, the HP approximation becomes exact 
and corresponding energy eigenvalues are double degenerate. This degeneracy is lifted for 
any finite number N of atoms by quantum tunneling between the wells of V±, a non-locality 
which cannot be captured by the HP approach. Moreover, the Hamiltonians derived within the 
HP approach are only suited to study the low energy sector of superradiant and normal phase, 
respectively. They definitely fail to describe the non-equilibrium dynamics when the system 
is initialized far from equilibrium, e.g. in a highly excited spin state, \n = 0,M = +N/2). 

To go beyond these limitations, an alternative formulation of the Dicke Hamiltonian 
[19] may offer some insight. For this purpose, one introduces collective operators a'.a with 
a = a + 2A/ (c&ftN) J x so that 


11 1 Dicke — ( 0 ( 1(1 


4A 2 

Nai 


(j*) 


2 + eft 


( 6 ) 


In a / eigenbasis, Dicke physics can then be visualized as governed by a set of shifted 
parabolas with nearest neighbor coupling described by the eft term (see figure 1). Starting 
for instance from the highly excited ft eigenstate \n = 0,M = +N/2), the corresponding 
quantum dynamics includes Landau-Zener-type of transitions through a multitude of avoided 
level crossings. Details of the corresponding time evolution will be discussed elsewhere. 


2.3. Non-equilibrium dynamics of the Dicke model 

The above discussion brings us to the question concerning to what extent the two facets of 
superradiant physics can appear within a single experimental setup. Here we argue that in 
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a minimal setting it suffices to add dissipation in the field mode and then to consider the 
relaxation dynamics of the compound from a maximally excited atomic state. 

In the regime of very strong damping, where the level broadening in the cavity mode by 
far exceeds the coupling energy between atoms and cavity, cavity excitations decay quickly 
and the emission process resembles that of a light burst into free space. The density of states 
of the open cavity leads to slight modifications in the transition rates between different spin 
states as compared to the emission into free space though. The corresponding dynamics of 
the collective spin after elimination of the cavity degree of freedom has been studied in depth 
previously [20], However, strong damping implies that the thermodynamics of the composite 
system may be strongly modified by the system-bath interaction so that no conclusions can 
be drawn from results found for the isolated system (e.g., the existence of a phase transition). 
In the opposite regime of vanishing dissipation, the dynamics of the Dicke model has been 
studied in the past with the focus on collapse-revival features, regularity, and quantum chaos 
[21]. This may give insight into collective properties of the compound but does not allow to 
access stationary state features for longer times. 

In the following, we thus concentrate on finite but weak damping, where the dominant 
dissipative mechanism is due to a large but finite quality factor of the cavity. In actual 
realizations of the Dicke model, where trapped ensembles of cold atoms are brought into close 
proximity with a superconducting cavity, competing dissipation mechanisms like spontaneous 
emission, trap losses, processes involving additional states of a single atom, or dark states of 
the atomic ensemble are expected to be of minor relevance only, although they may have 
an impact in specific ranges of parameter space [14]. For circuit QED realizations with 
artificial atoms other noise sources must be considered, with decoherence most prominent 
in the Cooper pair boxes. We leave this issue for future work and, in the spirit of a minimal 
setting, model the decay of electromagnetic excitations of the cavity in a standard manner by 
adding to the Liouville-von Neumann equation for the density operator of the bare composite 
system a Lindblad-type dissipator. Hence, we write for the time evolution of the reduced 
density operator of the cavity-atom aggregate 

^P = {[^Dicke,p]+^diss(p) (7) 

with the dissipator 

-^diss(p) = K"(« + 1) (2apa^ — a'ap — pa'a) + Kh(2a f pd — da'p — paa ( ). (8) 

Here, K is the damping rate and h is the thermal Bose occupation factor for an environmental 
photon with cavity energy ft). Note, that we employ the Lindblad dissipators derived for 
radiation damping of the cavity without coupling to the spin. This captures sufficiently 
accurately the non-equilibrium dynamics from the initial transient period up to long times 
where final equilibration sets in, and thus the regime where we expect the Dicke Hamiltonian 
to be realized in actual mesoscopic systems. Of course, for weak atom-field coupling, this 
modeling of dissipation is correct on all time scales. 
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In this section the dissipative dynamics of the field mode or, equivalently of the cavity 
occupation according to (7), is analyzed in detail on various time scales. The goal is to reveal 
the relation between and to characterize the nature of atom-field correlations during the initial 
and the final stage of the relaxation process. 

3.1. Emission burst and relaxation 

The situation is considered where initially the compound is prepared far from its equilibrium 
state with all atoms in the excited state and the cavity in its ground state. We first concentrate 
on the zero temperature case and later on comment on finite temperature effects. As expected, 
the atom-cavity system starts to develop a coherent flow of excitations associated with an 
oscillatory pattern in typical observables such as the mean cavity occupation (h) (see figure 2) 
and spin expectation values {J x / y / z ) (not shown). Due to the leaky cavity these features are 
damped out so that the system finally relaxes to a stationary state. In figure 2(a) we contrast 
the dynamics of (h) for couplings A below and above the critical coupling A c . Results for the 
RWA model are included as well. In the regime below A c , both the full Dicke model and its 
RWA version basically coincide and predict a decay towards (h) = 0. This is no longer the case 
in the regime above the critical coupling. There, the full model approaches asymptotically a 
finite cavity occupation in agreement with thermodynamic calculations (cf. [17, 18]) and (5)) 
which cannot be captured by the RWA dynamics. 

Insight into this failure of the RWA model can be obtained from the master-equation (7). 
One easily finds in the stationary limit ^ (h) l^oo = 0 = ^ (/ z ) |r->oo that 

(n) = n + i-({&?-) ~ > ( 9 ) 

where the latter counter-rotating terms are absent if the RWA is employed. Hence, within 
RWA the Lindblad-damping term used in (8) always drives the system into a the thermal state 
of the bare cavity independent of the cavity-atom coupling. Consequently, no signatures of the 
phase transition can be found in the dynamics of the cavity occupation for the RWA model. 
In contrast, in the full model (( aJ~) — (aU + )) plays the role of an order parameter which is 
zero below and takes finite values above the phase transition. Dissipation guarantees that the 
compound relaxes to the correct correlated thermal state which approaches with increasing N 
and for T = 0 the result known from the HP treatment. Of course, for any finite N the phase 
transition (for example (h) t ^ooVs. A) is smeared out. 

This crucial importance of the interplay of non-RWA terms and the particular damping 
mechanism for the relaxation towards the correct ground state has not been fully appreciated 
in previous studies of dissipative Dicke models [14, 22], Note though that in [14] the RWA 
together with dissipation induced by atomic losses also yields signatures of a phase transition 
in the dynamics. 

The benefit of the RWA is that the set of equations of motions for (h) and (J z ) close 
in the stationary limit. This is no longer the case for the full model. Analytical progress 
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Figure 2. (a) Dynamics of the mean cavity occupation (n) for coupling strength above/below 
the critical coupling and with/without the rotating wave approximation of the coupling term. 
Above the critical coupling and without RWA the system equilibrates to steady state with 
finite mean cavity occupation (for comparison the HP result is shown). Employing RWA 
(green dashed and blue dotted lines) the system does not equilibrate to its ground state but 
reaches a steady state, where the cavity is empty, (parameters: N = 10, k/co = 1/20, n = 0 
and A = 0.2 A c - normal, A = 2.4 A c - superradiant) 

(b) Occupation of eigenstates (arranged in order of ascending energy) for the initial state 
(\n = 0, M = +/)), an intermediate ’’burst” time (when (n) is maximal), and in the long¬ 
time limit. At the burst many eigenstates are involved. The steady state occupation highlights 
the potential structure above/below the phase transition (see subsection 3.1 in the main text), 
(parameters as above, but n = 0.1 for bottom right panel) 
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can be made, however, for the short-time dynamics starting from the initial state of fully 
excited spins. One finds for the initial buildup of cavity excitations in the weak damping limit 
(h) = A 2 t 2 (cf. figure 3) independent of the number of atoms N and valid for both the full 
model and for the RWA. 

To elucidate further details of the dynamics and particularly make progress in 
understanding the nature of the first emission peak in figure 2(a), it is instructive to analyze 
how different eigenstates |v) of the Dicke Hamiltonian participate during different stages of 
the dissipative time evolution. In figure 2(b) probabilities P v to find the system in eigenstate 
| v) at three different stages of the dynamics are displayed: at t = 0 (left), at the ’’burst time”, 
when the first peak in (h) occurs (middle), and for long times (right) for couplings below and 
above the critical coupling (RWA data are included as well $). 

While in the normal regime the low energy sector is dominantly occupied throughout 
the time evolution, this changes drastically for A > A c , particularly for times around the first 
emission peak. There, a large number of eigenstates participate substantially. Asymptotically, 
the probabilities P v show a near-degeneracy of pairs of eigenstates which is reminiscent of 
the exact degeneracy in the thermodynamic limit according to the HP treatment. As discussed 
above, finite size effects induce correction terms to the HP Hamiltonians H ± and lift their 
degeneracy. How many pairs of near-degenerate eigenstates occur, is a measure of the height 
of the barrier between the potential minima at (±\/2a, =F-y/2j8). To illustrate the appearance 
of these near-degenerate states, in figure 2(b) data are shown at a slightly elevated temperature 
h = 0.1. Note that one has also finite occupations in some higher lying eigenstates due to the 
small but finite value of the damping parameter k. Apparently, the data for the RWA model 
strongly deviate from these findings for A > A c . In particular, in the stationary limit the 
system is not relaxed to the lowest energy eigenstate but to a state close to the ground state of 
the uncoupled system. 

3.2. Characteristics of the emission burst: Is it superradiance? 

In this subsection we will argue that the first pronounced peak in the cavity occupation of 
the damped Dicke model carries signatures of the superradiant free-space burst and reflects 
its characteristic features. These are the crossover in the time-dependence of the emission 
rate from an exponential decay for a single atom to a peak, where for larger N the height 
of the peak grows as N 2 , its temporal width decreases as N~ l and the integrated emission is 
proportional to the total initial excitation energy and, hence, to N. 

These scaling properties are modified when one considers a fixed value of the parameter 
A upon increasing N, the typical scenario for the analysis of the phase transition in the Dicke 
model. Then, the rescaled coupling constant between a single atomic excitation and the e-m 
field is g = A/ y/N so that the emission rate due to a single excitation is reduced g 2 = A 2 /N 
for a larger number of atoms. Accordingly, the characteristics of the light emission are as 
follows (see figure 3): (i) For N independent atoms the radiation follows an exponential decay 
from a fixed (/V-independent) value over a timescale growing linearly with N; (ii) For a super- 

$ Far below the critical coupling the probability distribution with and without RWA is obviously very similar. 
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radiant burst of N cooperating atoms a peak with height N and constant width °< N° emerges. 


We can now test these signatures for the mesoscopic dissipative Dicke model (figure 3). 
Indeed, for A > A c we find an increase of the height of the first peak with growing N, while 
its width remains basically constant. Quantitatively, for the ^/-dependence of peak height and 
width (defined here as the left-sided half width at half maximum value) the expected scaling 
with N and N°, respectively, is obtained (not shown). While for the results discussed so far 
we considered the resonant case, where the cavity energy equals the energy needed for atomic 
excitations, we numerically checked also the off-resonant case. It turns out that the discussed 
dynamical features are not very sensitive to weak de-tunings. This is in stark contrast to 
the purely coherent dynamics, where complex collapse-revival dynamics with domains of 
regularity and quantum chaos occur [21]. We enter this dynamical range for extremely weak 
dissipation as well. 

In conclusion, our findings strongly suggest that the physics of the free-space 
superradiant burst can indeed be recovered in the damped dynamics of the Dicke model 
for suitable parameters. In an actual experimental realization with a cold-atom cloud 
coupled to a superconducting stripline-cavity along the lines of [10], the atom number N 
is orders of magnitude larger than values which can be realized in full numerical simulations. 
Nevertheless, the scaling characteristics of the peak will allow for a wide range of damping 
rates, where friction is strong enough to suppress coherent collapse-revival dynamics but weak 
enough to achieve prominent radiation bursts with peak cavity occupations 1. 

4. Monitoring collective dynamics 

4.1. Dicke dynamics in phase space 

Our discussion of the dynamics of the dissipative Dicke model has so far been focused on 
a single observable, the mean cavity occupation (h), which is certainly not the only way, to 
gain experimental insight into the dynamics. A more complete characterization of the cavity 
dynamics is offered for example by the (reduced) Wigner density of the cavity degree of 
freedom, which is routinely measured for stripline resonator systems (cf. [23] for a recent 
example). 

In particular, in the thermodynamic limit N —> °° of a nearly classical spin, we may expect 
the Wigner density, which has a simple classical interpretation, to be exceptionally suitable. 
One question of interest is to distinguish sub- and superradiant coupling without waiting for 
the relaxation of the system to its stationary state. That way, such a distinction would be 
immune to the danger that the long-time limit is influenced or covered by extraneous terms, 
not present in the Dicke Hamiltonian (e.g., non-homogenous coupling [12, 13]) or additional 
or different damping mechanism (e.g., atom decay [14]). In fact, we find that the Wigner 
density dynamics shows clear signatures, reflecting the contrast between a unique mean field 
solution for the ground state below the critical coupling and two equivalent solutions above, 
already for short times and a small number of atoms. 
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Figure 3. Dynamics of the mean cavity occupation (n) of the superradiant damped Dicke 
model (A = 1.2Ac, k/o) = 1/40) for increasing number of atoms, N = 5,10, 15,20. 

The insets contrast the emission of light into free space from N independent atoms (normal) 
and collective emission (superradiant) where the single-atom to field coupling is rescaled 
g ~ 1 / fN as in the Dicke case. The first peak in (n) has the same scaling characteristic 
as the superradiant burst into free space: the peak height increases linearly with N, the width 
is constant. 


In figure 4(a) we show the reduced Wigner density of the cavity (in dimensionless 
coordinates ( X a ,P a )) well above the critical coupling at various times. The Wigner function 
of the initial state at t = 0 is a Gaussian peak centered at (X a = 0, P a = 0). It starts to spread 
out in two arms in the positive and negative X a direction and evolves into the structure shown 
in figure 4(a) ii) at the time of the burst, t = /burst- A rich fine structure of complex interference 
patterns remains for some time, while substantial weight is assembling in two peaks (cf. 
(a) iii at / = 2/burst)- These peaks rotate around and finally slowly relax towards two final 
positions in phase space (cf. (a) iv), which correspond to the two equivalent solutions of the 
HP approximation in the N —>■ °° limit. 

The spin-dynamics can similarly be visualized by another quasi probability-distribution, 
the Husimi function, relying on the concept of spin-coherent states. We show the final state’s 
reduced spin-Husimi function in the phase space spanned by angular variables (0,0), which 
indicate the ’’direction” of the spin vector, in figure 4(b) for the same superradiant coupling 
as in (a). As for the cavity degree of freedom, we find a two peak structure for spin values 
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Figure 4. Dynamics in phase space. 

(a) Reduced Wigner density of the cavity above the critical coupling (A = 1.2Ac): Starting 
from a ground state, i, the spin system excites the cavity. Occupation oscillates around, ii at 
t = tburst and iii at t = 2/b urst , and finally relaxes to two minima positions iv. 

(b) shows the Husimi function of the spin degree of freedom in the stationary state, which 
shows a similar relaxation towards two final peaks. The spin is then directed towards two 
opposite points on the southern hemisphere (cf. the equator indicated by the white circle, 
while the south pole (6 = n—6 = 0, corresponding to \M = — J )) is in the center of the plot. 

(c) shows, that the spin-reduced M x = —J component of the cavity Wigner density relaxes to 
one of the stationary state positions. 

Below the critical coupling (A = 0.2A C ) the cavity Wigner density starting from the ground 
state, cf. (a) i, undergoes a breathing-type oscillation, (d) i at t = fb U rst> and finally relaxes to a 
final peak close to the ground state, (d) ii. (other parameters: N = 6, k/w = 1/10, n = 0) 


corresponding to the two equivalent solutions of the HP result, cf. (5). 

For subcritical coupling the initial central Gaussian peak broadens [see figure 4(d) i for 
X = 0.2 A c and t = /burst] and re-sharpens in a breathing-type oscillation, which is damped out 
to relax to a final Gaussian peak centered around the unique minimum (0,0) [see figure 4(d) 

ii]. 

The long-time behavior observable in figure 4(a) iii, iv and (b) can be intuitively 
understood in the HP picture as damped dynamics well within the two equivalent parabolic 
potentials in cavity and spin variables. For the short-time behavior, however, the HP 
approximation is not applicable. To explain the observed splitting of the central peak into 
two arms [figure 4(a) ii], one can employ the representation of the Hamiltonian in a J x basis 
yielding the shifted coupled parabolas of equation (6). Starting in the initial state, which is not 
a J x eigenstate, each J x component (which - as a cavity ground state - is a centered Gaussian 
in x-space) evolves in its correspondingly shifted parabola and starts to move into positive 
or negative x-direction. According to this picture, we expect, e.g., the M x = ±J component 
to start to move into positive/negative x-direction. Indeed, this behavior is observed in the 
corresponding spin-projected Wigner density of the cavity in figure 4(c) [same coupling and 
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time as in (a) ii]. This argument has so far neglected the coupling between the motion in 
the various parabolas induced by the J z term in (6), which causes transitions and interference 
structures, which become prominent, when the various component-wavepackets meet at a 
crossing of the parabolas. 


4.2. Indicators of sub/superradiant coupling 


The concept of a phase transition applies to a system in the thermodynamic limit (N —> °°), 
which after a long time it —» °°) is driven to equilibrium by a vanishingly weak coupling to an 
environment (k —> 0). For that case, the mean cavity occupation, (h)/J, is an order parameter 
for the Dicke model, i.e., an observation of its value clearly indicates whether the coupling 
is sub- or superradiant, A ^ A c . In our case, however, we are interested in an observable 
distinguishing between the weak and strong coupling case for a mesoscopic (N > 10 — 100) 
Dicke aggregate with finite coupling to the environment on short to intermediate time scales. 
The latter requirement is due to the presence of additional physical processes in any actual 
experimental realization, which our description of the dynamics of the damped Dicke problem 
does not take properly into account as discussed in the Introduction. 

It is, therefore, advisable to search for alternative observables, which may in this 
scenario be better indicators of sub- or superradiant coupling than the cavity occupation. The 
investigation of the dynamics of the Wigner density has revealed one distinctive difference 
between the sub- and superradiant case already apparent after short times: namely, a 
substantial fragmentation of the Wigner density into two parts, oscillating and finally relaxing 
to the two equivalent minima of the superradiant case. This is, for instance, reflected in the 
curvature of the Wigner density (in X a -direction) 




=0 =P a 


-- n nP 2 a) 

71 


at the center of the cavity phase space X a = 0 = P a , where the second expression containing 
the parity operator for the cavity degree of freedom, tl n , follows from the definition of the 
Wigner density. 

Note, that this observable, which can for example be gained from a subset of 
measurements of the full Wigner density, is but one example of a higher order correlation 
function of cavity operators. Other combinations, also mixed cavity-spin correlators may 
be of similar or in other parameter regimes even of superior use for deducing the coupling 
strength from observations of the short-time dynamics. In figure 5, we demonstrate that the 
sign of the new observable can indeed give a clear indication of sub- or superradiant coupling 
strength after a short observation time of the order of the inverse damping rate, e.g., Kt < 3 
(lower row). In contrast, the dynamics of the order parameter of the thermodynamic case, 
(h) /J, up to this time only does not easily allow to gauge, how close to relaxation the system 
is (upper left panel of figure 5); even less does it show any distinguishing features around the 
transition from sub- to superradiant coupling (symbols in upper right panel). 
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Figure 5. Observables as indicator of normal or superradiant phase. 

The order parameter, (h)/J, of the thermodynamic phase transition does not indicate clearly, 
if the coupling k is above or below the critical value A c , when the dynamics of the damped 
(mesoscopic, N = 6) Dicke model is observed on short time-scales only (upper row). 

Other observables, like the curvature of the (cavity-)Wigner density, (TL n P 2 ), which quickly 
goes to positive (negative) values for coupling values below (above) the critical value, can give 
clearer signatures (lower row). 

The left column shows the dynamics of the two observables for various values of the coupling 
(+ : A/A c = 0.8, ♦ : A/A c = 0.96, x : A/A c = 1.12, ■ : A/A c = 1.28,) for times on the order 
of the inverse damping rate, Kt ~ 1. 

In the right column the value of the observables at an intermediate measurement time, Kt = 3 
(symbols), is contrasted to its stationary limit with (k/co = 1/10 - magenta), and without 
damping (black solid line). The order parameter observable, (h)/J (upper right panel), 
while distinctive in thermodynamic limit (black dashes), is less suited as an indicator of 
super/subradiant coupling, than other observables, like the Wigner curvature (lower right 
panel), which shows a distinctive sign change at (approximately) the critical coupling strength. 


5. Conclusion 

Investigating the dissipative Dicke dynamics in a system where a mesoscopic cold atom 
cloud is coupled to a superconducting cavity, we found that two facets of superradiance 
physics appear. In the long term, dissipation drives the relaxation of the system towards 
its equilibrium state, which shows the profoundly different properties of a normal and a 
superradiant phase below and above the critical coupling strength. On short time scales, a 
peak in the cavity occupations was shown to exhibit the characteristic scaling properties of 
the free-space superradiant radiation burst. The radiation burst is further characterized by the 
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participation of many eigenstates in the dynamical evolution for which non-RWA terms are of 
crucial importance. 

To overcome issues problematic for observing the long-time limit of full relaxation, we 
investigated the phase-space dynamics of cavity and the atomic degree of freedom, where 
the reduced (quasi-)probability distributions were found to reflect the appearance of two 
degenerate minima in the superradiant regime. This allows to identify suitable obervables and 
features in their transient dynamics, which allow a clear distinction of sub- and superrradiant 
coupling on short time scales. 
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